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Orthogonality catastrophe in fermionic systems is well known: in the thermodynamic limit, the 
overlap between the ground state wavefunctions with and without a single local scattering potential 
approaches zero algebraically as a function of the particle number N. Here we examine the analo- 
gous problem for bosonic systems. In the homogeneous case, we find that ideal bosons display an 
orthogonality stronger than algebraic: the wavefunction overlap behaves as exp[— XN 1 ^ 3 ] in three 
dimensions and as exp[— XN/ In 2 N] in two dimensions. With interactions, the overlap becomes 
finite but is still (stretched-)exponentially small for weak interactions. We also consider the cases 
with a harmonic trap, reaching similar (though not identical) conclusions. Finally, we comment on 
the implications of our results for spectroscopic experiments and for (de)coherence phenomena. 



In the area of correlated electrons, there is a long his- 
tory of studying quantum impurity problems such as the 
Kondo effect A broad range of phenomena arise de- 
pending on the type of impurity, the nature of the bulk 
electron system, and the way they are coupled. The con- 
ceptual basis for this richness is provided by the orthogo- 
nality catastrophe Q . It deals with the effect of a single 
scattering potential on the many-body states of an ideal 
iV-fermion system. The many-body ground state wave- 
function is a Slater determinant of the N lowest single 
particle wavefunctions. Since the impurity potential af- 
fects each and every one of the single-particle states, its 
effect on the many-body state is significantly amplified: 
in the thermodynamic limit the ground state wavefunc- 
tion in the presence of the impurity potential (|gs' >) 
is orthogonal to its counterpart in the absence of the 
impurity potential (|gs >). For large but finite N, the 
wavefunction overlap, S = | < gs|gs' > |, depends on N 
in an algebraic form: 

Sfcrmion ~ AT . (1) 

The exponent a = 8 2 /dir 2 when the scattering contains 
only an s-wave component; here, 8 is the s— wave scat- 
tering phase shift of the single-particle state at the Fermi 
energy, and d the dimensionality The initial work of 
Anderson provided an upper bound. An exact solution 
was later achieved in the work of Nozieres and collabo- 
rators Q. For an illuminating discussion, including the 
connection with the Friedel sum rule, see Ref. Q. 

Cold atoms provide a setting to engineer many-boson 
systems with a variety of quantum phases |f| . It is natu- 
ral to ask what happens when local defects are introduced 
into these systems. Here we consider the orthogonality 
effect in a Bose-Einstcin condensate (BEC). For ideal ho- 
mogeneous bosons, the problem will be solved exactly in 
a very simple way. For the cases with interactions and/or 
confining background potentials, we will carry out anal- 
yses pcrturbativcly in the impurity potential. 

Ideal bosons in a uniform background: In this case, 
the effect of an impurity can be determined exactly. The 



ground state condensate wavefunction is 

|gs >= ^o(xi)<ySo(x 2 )---<A)(xjv), (2) 

where ipo is the single-particle ground state wavefunction. 
|gs' > has the same form, with <y9 being replaced by 
ip' , the single-particle ground state wavefunction in the 
presence of the impurity. The overlap, S, is then simply 

S = s N (3) 
s = < tp \tp' > ( 4 ) 

Consider first the case of three dimensions. We choose 
a spherical box of radius R, and a fixed boundary condi- 
tion ip(r = R) = 0. The impurity is placed at the center: 

AH = VQ{a-r) (5) 

where a is its size. Without a loss of generality, the 
potential will be taken as repulsive, we define ip as 
the radial part of the single-particle wavefunction multi- 
plied by r. In the region < r < a, ipi — Asin(fc/r), 
with fc/ 2 = 2™(e - V). In the region a < r < R, 
ipjj = Bsin(knr + 8), with fc// 2 = ^pre, where 8 is the 
phase shift. The fixed boundary condition at r = R 
yields kjj = ^j^-- The continuity of ip and its derivative 
at r = a gives rise to 

tan(fcja) tan(fc//a + 8) 

1 = 1 ( 6 ) 

fc/ fc// 

The phase shift for the single-particle ground state in 
the thermodynamic limit (i?/a — > oo) is found to be 

8 3d =~f(V) (7) 

where f(V) = tt(1 - iSHM£i ) ? wit h V = (^^) 1 / 2 . The 
a/R dependence reflects the fact that the single particle 
state sits at the edge of the band, where the density of 
states is proportional to ^ft oc a/R. The numerical result 
[Fig- El shows the validity of this form already for R/a of 
the order 10 3 . 
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From the forms of the wavefunctions in the presence 
and in the absence of the impurity, we determine the 
overlap of the single-particle wavefunctions 



(8) 



where we have used the fact that 63d is always small in 
the thermodynamic limit. Combining Eqs. 1318171) . we 

find a stretched-exponential form for the overlap of the 
many-body ground states: 



S 3d = exphA(^)iV 1/3 ] 



(9) 



( Aung ' 
k 3 , 



with uq being 



where \(V) = + -^)a 2 f 2 {V) 
the particle density. 

Consider next the case of two dimensions - a disk of 
radius R. The ground state wavefunction ip is now a 
linear combination of the zeroth order Bessel and Neu- 
mann functions: ipi = AJo(kir) and ipn = B[Jo(kur) + 
tan((5)TVo(fc//r)]. The fixed boundary condition at r = B 

*-+5 



now gives kji 
its derivative leads to 



R 



The continuity at r = a of ip and 



kjJi (kja) kjj [J 1 (kjja) + tan(<5)TVi (kjja)] 



J (kia) 



Jo(kna) + tan(<5)A r o(fc//a) 



(10) 



Using the limiting forms of the Bessel and Neumann func- 
tions appropriate for a/R — > 0, we find the following 
phase shift of the single-particle ground state 



1 



J2d 



2 ]n(R/a) 



(11) 



Note that it is independent of the potential strength. The 
finite value of the density of states at the band edge would 
have implied a phase shift that is independent of R/a. 
The exact result, on the other hand, contains a logarith- 
mic correction factor. In Fig. 2, we plot the numerical 
solution which confirms the logarithmic factor. 

Taking these results together we find the overlap of the 
many-body ground states in two dimensions: 



5 2d = cxp[-aTV/ In 2 TV] 



(12) 



where a = 4r 



' f J a (3j-x)B(xy 



with 



B(x) = N (^-x) - Ji(^fx)x and / /(*) = ft f(x)xdx. 

So far, the results are exact. In order to generalize to 
situations in which exact solutions are not readily avail- 
able (see below), we now turn to a perturbative treat- 
ment of the impurity potential. To the second order in a 
linked cluster expansion, the overlap between the many- 
body ground state wavefunctions is jf|: 



S 
C 



cxp[— C] 

1 I < «|A#|gs > I 2 



(E n - E gs ) 2 



(13) 



The impurity scattering can be rewritten as AH = 
Va'{~x = 0)a(x = 0), where a'(x) and a(x) are the 
boson field operators. The off-diagonal long-range or- 
der of a BEC implies that we can set ao and at to 



be VA^. The impurity Hamiltonian becomes AH = 

AHi + V/L d {N + E p ^p'«iv). where the P art of 
concern to us is 



(14) 



Here, L is the linear dimension of the system. Combining 
Eqs. (|13I14|) . and using the dispersion relation e(p) = 
p 2 /2m, we have 



C = 



V 2 N (2m) 2 
T17 {2nh) d 



d d p p~ 



Htt/L 



(15) 



In three dimensions, C has a divergence of the form 
L cx iV 1 / 3 , which is the same as the exact result [Eq. 
(For the spherical box geometry considered earlier, the 
perturbative result is found to fully agree with what 
comes out of Eq. (JSJ when X(V) is expanded to order 
V 2 .) In two dimensions, the divergence becomes L 2 cx TV. 
Compared with the exact result [Eq. ^1 > the perturba- 
tive result recovers the TV factor in the exponential but 
misses the multiplicative logarithmic correction. 
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The pcrturbative treatment also provides the physical 
picture for our results. To see this, we rewrite the expres- 
sion for C in terms of the dissipation spectral function, 
Jm%Q as follows, 

^Imx^H (i6) 



ImXoV) 



< n|Ai?|gs > 



xS(lu - E n + E gB ) 



(17) 



where E\ is the energy of the first excited state with a 
non-zero < n|Ai?|gs >. The form of the impurity poten- 
tial, Eq. (|14fl implies that the dissipation spectral func- 
tion is simply proportional to the single-particle density 
of states: 



Imxo 1 (w>0)aw2 1 . 



(18) 



The exponent (| — 1) reflects the quadratic nature of the 
dispersion of the low- lying excitations. It is less than 1 
for both three and two dimensions; in the terminology 
adopted in the dissipated two-level system literature [7J, 
both arc sub-ohmic || . This abundance of low-lying ex- 
citations is responsible for the strong orthogonality in a 
BEC. 

Weakly interacting bosons in a uniform back- 
ground: Unlike for fermions, even weak interaction is 
a relevant perturbation for bosons. We use the standard 
Bogoliubov transformation 0, 0| , 



1 



L p bl 



(19) 



Under this transformation, the quadratic part of the 
Hamiltonian becomes 



p#0 



(20) 



For p 3> 2mu (where u = ^JUno/m is the sound 
velocity, with U being the effective contact interaction 
amplitude), L p nearly vanishes and we recover the non- 
interacting limit, including Eq. (|14fl and£(p) « e(p). For 
p <C 2mu, on the other hand, L p f=a —1 +p/mu, and it 
follows that 



£(p) 



up 

p#0 v 



where the prime denotes that the summation is up to 
about 2mu. Eq. (|15|) is then replaced by 



C 



V 2 N 1 ,„ , 2 

^l^(2^F (2m) x 



(2mu) 



I p2mu poo 

— / d d P p- 1 + / d d P P - 

%u ) JUtt/L Jlrnu 



(22) 



It is convergent in both three and two dimensions, imply- 
ing a finite wavefunction overlap. The resulting overlap 
in the thermodynamic limit depends on the interaction 
U in the following (stretched-)exponential forms: 



Su,3d = cxp[-X'/VU} 
Su,2d = exp[-a'/U] 



(23) 



where A' 



,:i/2 



47T 2 fi 3 



and a' 



3V* m 



This conclusion can also be seen through the form 
of the dissipative-bath spectral function, which, at low- 
energies, now takes the form [cf. Eqs. i|17l21|) ] 



Imx > 0) (x lu° 



(24) 



Its super-ohmic nature in two and three dimensions im- 
plies that C [cf. Eq. (|16|) ] is infrared convergent 0. 

Ideal bosons in a harmonic confining potential: Con- 
sider an isotropic harmonic trap with frequency Qq. The 
thermodynamic limit is defined by keeping 



const. 



(25) 



as both N and l/f2o go to infinity [llj]. (For instance, it 
ensures a finite energy per particle for ideal fermions.) 

We consider an impurity located at the center of 
the trap: AH = Vd{x = 0)a(x = 0). Using a(x) = 
^ n <^ n (x)a n , where n = (rii,- • ■,rid), and the single- 
particle eigenfunctions 



Vn(x) 



2h 



d/4 



ntrK!)- 
n d 



i / mQo 



2h 

mSlor 2 /2h (26) 



mflo dx 

we write the linear part of the impurity Hamiltonian as 



AHi = V^Nn a d J2^2 n (al 

n^O 



+ a 2n ). (27) 



Here, we have defined 



_^ n (0Vo(0) (-i)£SU 

^2n = 



n 



(2n, - 1)! 
(2rn)!! 



From Eq. I|13|) . we have 
V 2 



c 



1 



2 m "\m y^ o {2nr 



E 



A 2 



(28) 



(29) 



where n = y ■_ 1 rij. At large re, we find that 

J2 ni ,-,n d S ( n ~ J2t=i n i) A 2n is proportional to v/n in 
three dimensions and approaches a constant in two di- 
mensions 0], so the summation over re in Eq. (|29|1 is 
convergent at ultraviolet. C is then proportional to f2^~ 2 , 
which, using the thermodynamic limit, is equivalent to 
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N 2 / d . The overlap between the ground state wavefunc- 
tions is then 

SWsd oc exp[-A /io 7V 2/3 ] 

Sho,2d oc exp[-a;, 7V] (30) 

where X ho = ^C,{\)(m/irh 2 f{Nh z 9, )^ and a ho = 

V 2 (m/h 2 ) 2 /AS. 

Weakly-interacting bosons in a harmonic confining po- 
tential: To understand the effect of interactions, we 
first note on one important consequence of the thermo- 
dynamic limit, Eq. I|25|) . In the non-interacting case, 
it is seen, from the ground state wavefunction ^o(x) = 
(m,n /irh) d/i e - mQ ° r2 / 2h , that the central density of the 
condensate, n(x = 0) = N\(Pq(x — 0)| 2 , is of the order 
of \/N in the thermodynamic limit. 

In the interacting case, on the other hand, n(x = 0) 
is well-known to be of order unity in the thermodynamic 
limit [ill ]. In this limit, the interaction term dominates 
over the kinetic term |l3| . It follows from the Gross- 
Pitaevskii equation that n(x = 0) = [ijU where the 
chemical potential [i is finite in the thermodynamic limit. 
This implies that <^o( x = 0) in the interacting case is a 
factor of iV -1 / 4 smaller than its counterpart in the non- 
interacting case. The wavefunctions of the low-lying ex- 
cited states should contain a similar reduction factor. On 
the other hand, the energies of the collective modes re- 
main linear in Q . It follows [cf. Eqs. I|13I27I28|) ] that 
interactions weaken the orthogonality effect, as in the 
uniform case. The precise form of the overlap depends on 
the details of the excited-state wavefunctions, which can 
be determined from the Bogoliubov-de Gennes equations 
in a harmonic potential; this will be discussed elsewhere. 

Experimental implications: In addition to the theo- 
retical significance, the orthogonality effect may also be 
directly probed in experiments. One way is to perform 
the analog of the x-ray edge measurement in metals 0,0- 
Consider a condensate co-existing with a separate species 
of atoms that are considerably more dilute and are lo- 
calized (by a deep optical potential well that only these 
atoms see). The photo-absorption or luminescence cor- 
responding to a transition between two levels of this sec- 
ond species of atoms would ordinarily be a sharp delta 
function. (In practice, the spectral width of a hypcrfinc 
transition for atoms in a BEC can be as narrow as ^100 
Hz or even ~10 Hz 0]-) However, the two atomic levels 
will provide different scattering potentials for the atoms 
of the condensate. The weight of the delta function - 
which is precisely the overlap of the condensate wave- 
functions corresponding to these two different potentials 
- would then have to vanish in the thermodynamic limit 
due to the orthogonality catastrophe. A finite spectral 
weight can arise only when the condensate atoms go to 
the excited states under the new potential, where the ex- 
citation energy serves as a cutoff for the infrared diver- 
gence. This results in a one-sided spectrum which can 



be divergent or vanishing at the edge depending on the 
degree of orthogonality for the low-lying excited states. 

The orthogonality may also be manifested in the time 
evolution of a condensate after a sudden introduction of 
a local potential. The orthogonality makes it rather hard 
for the system to evolve into the new ground state. In 
other words, the density distribution will tend to keep 
its initial profile; the impurity is "hardly visible" to the 
condensate. 

Yet another implication is on the coherence and de- 
coherence phenomena. Consider, for instance, localized 
atoms in a BEC. An effective Kondo problem arises when 
both the localized atoms as well as the itinerant atoms 
of the condensate contain real or pseudo- spin degrees of 
freedom 0, Strong orthogonality makes the spin 

flip process entirely incoherent. 

To summarize, we have studied the orthogonality effect 
in Bose-Einstein condensates. For ideal bosons, the over- 
lap of the ground state wavefunctions when a single local 
scattering potential changes from one value to another 
vanishes in a strctched-exponential form in the thermo- 
dynamic limit. With interactions, the overlap becomes 
finite but is small for weak interactions; its dependence on 
the interaction strength is typically stretched-exponential 
as well. These effects can be probed using spectroscopic 
experiments in cold atoms, which can be tuned from be- 
ing essentially ideal to strongly interacting. The effects 
also have significant implications for the coherence and 
decoherence phenomena in bosonic systems. 
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